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THREE-DIMENSIONAL POTENTIAL AND CURRENT DISTRIBUTIONS IN A 
HALL GENERATOR WITH ASSUMED VELOCITY PROFILES 
by N. Stankiewicz and R. W. Palmer 
Lewis Research Center 

SUMMARY 

Three-dimensional potential and current distributions in a Faraday segmented mag- 
netohydrodynamics (MHD) generator operating in the Hall mode are computed. Constant 
conductivity and a Hall parameter of 1.0 is assumed. The electric fields and currents 
are assumed to be coperiodic with the electrode structure. The flow is assumed to be 
fully developed and a family of power-law velocity profiles, ranging from parabolic to 
turbulent, is used to show the effect of the fullness of the velocity profile. Calculation 
of the square of the current density shows that nonequilibrium heating is not likely to 
occur along the boundaries. This seems to discount the idea that the generator insu- 
lating walls are regions of high conductivity and therefore responsible for "boundary- 
layer shorting", unless the shorting is a surface phenomenon on the- insulating material. 


INTRODUCTION 

Two-dimensional potential and current distributions within a magnetohydrodynamic 
(MHD) generator channel have been investigated in many studies (see, e.g. , refs. 1 
to 4) . These two-dimensional results are adequate if the velocity profile is constant in 
the magnetic -field direction. Unfortunately, the viscous effects of the side walls cannot 
be ignored in a real MHD channel because of its finite width. 

It is suspected that the shorting noticed in real generators is caused by the velocity 
profiles in the channel. Near the insulator sidewalls, the smaller stream velocity 
causes a reduction in the Faraday component of the current (see fig. 1). This, in turn, 
reduces the f x S force against which the stream does work. A Hall potential is induced 
in the downstream electrodes as the result of this work. For a given Hall current, the 
Hall potential will therefore be less than the expected value without insulator sidewalls. 



Conversely, for a given Hall potential, the Hall current will be higher as a result of the 
sidewalls, giving the appearance of a short circuit. 

Because of the importance of the shorting problem in the operation of an actual gen- 
erator, it is necessary that the effect of the sidewalls be studied. This necessitates a 
three-dimensional analysis. 

The problem studied in this report is a constant conductivity, fully developed flow 
problem in which the velocity profile varies only transverse to the flow. It is shown in 
reference 5 that a three-dimensional analysis is not necessary if one is only interested 
in the overall duct characteristics (load lines); that is, by averaging over the third di- 
mension (parallel to the magnetic field) the problem can be reduced to a two-dimensional 
problem. This method, however, does not determine the local potential and current dis- 
tributions. It is the local current distribution that affects the local conductivity and 
efficiency. 

The purpose of this study is to solve the three-dimensional potential problem when 
the velocity profile varies across the duct. However, rather than solve the coupled mo- 
mentum equation for the actual velocity profile, a power -law variation is assumed. That 
is, the velocity varies as a power of the distance from a wall. This is valid because the 
interaction parameters between the plasma equations and the fluid equations are small. 
The essential features of a three-dimensional problem are retained with this assumption. 

The exponent in the velocity function is a parameter in this study . A value of 1 . 0 
for the exponent yields a parabolic velocity profile; and a value of 1/7 is a good approxi- 
mation to turbulent flow across a flat plate. For turbulent flow in a channel with a square 
cross section the power seems closer to 1/6 (ref. 6). Varying the velocity profile 
changes the effective size of the region in which shorting is expected. 


SYMBOLS 


A(K,L,M) 


a 

B 

b 

C 

D 

E 

e 


computational definition of normalized potential function used in finite dif- 
ference forms of differential equations 

dimension of duct in y direction 

magnetic field 

dimension of duct in x direction 

number of points into which a unit line is divided for computational use 
damping factor 
electric field 
electron charge 
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(x,y,z) 


Greek Symbols 

Hall parameter 
ratio h/a 
ratio h/b 

electron-ion collision frequency 
normalized coordinates, (£ | |B; t ]\ | u; £ | |B x u) 
conductivity 
potential function 
normalized potential function 

Math Symbol 

( ) average value of 


|3 

r 

^ei 

a 

<P 

* 


dimension of duct in z direction 
total current 

normalized total current density 
current density 

integer point location in normalized x direction 
Hall load factor 
constant (eq. (4)) 

integer point location in normalized y direction 
integer point location in normalized z direction 
electron mass 

exponential parameter in velocity distribution 
time 

normalized velocity 
velocity 

present value of normalized potential function 
coordinates, (x||B;y||u; z||Bxu) 
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Subscripts: 
i ion 

e electron 

F Faraday 

H Hall 

(x , y , z) ve ctor c omponent 

Superscripts: 

— vector 

unit vector 

i iteration number 


ANALYSIS 

Derivation of Equations 

The fundamental equations used in the analysis of an MHD generator are the con- 
tinuity of current 


div ] = 0 


( 1 ) 


and the Ohm’s law relating current and electric field through a vector conductivity. 
Equation (1) can be written in terms of a scalar conductivity as 

j* = ct(E + uXB) — f x B 

My. 
e ei 

In addition, the basic equation of electrodynamics curl E = 0 allows the defining of E 
with a scalar potential <p, 


( 2 ) 


E = -grad <p (3) 

The coordinate system and field orientation are shown in figure 1. The flow is perpen- 
dicular to the magnetic field. 
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Figure 1. - Coordinate system and channel orientation. 


The velocity vector u is chosen in the positive y direction and a constant magnetic 
field B in the positive x direction. The pressure, temperature, conductivity, and 
channel cross section are held constant. End effects of the channel are ignored. 

A power-law velocity profile of the form given in equation (4) is assumed. The pro- 
file is constant in the y direction. 


u = K o [x(b - x) ] 1/m [z(h - z) ] 1/m (4) 

Equation (4) yields a parabolic profile when m = 1 and a turbulent profile when m •= 7. 

Equation (4) can be normalized to the average velocity (u) by means of this defini- 
tion: 

<u> = ^2 / b / h dx dz [x(b - x)] 1/m [z(h - z)] 1/m (5) 

kh O£=(Kz=0 
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Eliminating K with equation (5) , equation (4) becomes 


u = (u)bh 


/7 h 

•'O • / 0 


[x(b -x)] 1/m [z(h - z)] 1/m 
[x(b - x)] 1/m [z(h - z)] 1/m dx dz 


(6) 


Taking the dot product, and cross product of equation (2) with B and solving for 
j* • B and j x B one can solve the identity B x (j x B) = j B 2 - B(j • B) to get 


] = 


i + p L 


E + u X B + /3Bu + B X E + B(E • B) 


B 


B‘ 


(7) 


where the condition U • B = 0 is used, and /3 = eB/M e ^ e . is the Hall parameter. 
Using equation (3), equation (7) in component form becomes 


„ d(fi 

J x = ~ a — 

x dx 


Jy = 


J + s 2 V 3y Hz J 
: ” Li£-3£. uB l 

! 2 V 3y dz / 


( 8 ) 


( 9 ) 


( 10 ) 


Using equations (8) to (10) and remembering that u = u(x, z), equation (1) becomes 


(1 


£j2) CTJ£_ 

dx 2 


+ d\ L+ dh L = 


3y 


Sz 


. B — 

9z 


(ID 


The solution of equation (11) with suitable boundary conditions will give the potential 
distribution within the MHD channel. 

The Hall current Ijj is defined as the total current that passes from one unit cell 
a x b x h to the next unit cell. The surface integral of the current density normal to the 
plane y = 3a/4 (the midplane of the electrode insulators) is taken as the Hall current. 


6 



Hence, 




(12) 


Except for the numerical difficulties of integrating over the rapidly changing fields 
near the electrodes, each plane y = constant along the electrode insulators (a/2) < y < a 
should give the same value as equation (12). This can be shown by integrating div j 
over a volume element bh Ay (i. e. , two adjacent planes). 

The Faraday current I F is defined as the total current emitted by one of the elec- 
trodes in a unit cell. This current must be equal to the current flowing in the external 
Faraday circuit of the generator. In order to avoid the numerical difficulties when inte- 
grating near the electrodes, the Faraday current is calculated at the midplane of the cell 
z = h/2. 


% " 


y ’ !) 


(13) 


Again it can be shown, by integrating div j over adjacent planes z = constant, that 
equation (13) should be equal for each plane. 


Boundary Conditions 

In the Hall mode of generator operation, opposite pairs of electrodes (opposite in the 
Faraday sense) are shorted, and the potential is allowed to build up in the Hall direc- 
tion. Referring to figure 1, the Faraday electrodes are located at z = 0 and z = h with 
y = a/2. (The lines located at z = 0 and z = h and at y = a are the leading edges of 
the next pair of electrodes.) 

The potential of the first pair of electrodes is taken to be zero, and this provides a 
Dirichlet condition for the problem; that is, 

cp = 0 at z = 0,h (14) 


0 <y 

2 
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The current density normal to the sidewall insulators located at x = o and x = b 
mush vanish. Using equation (8) yields the Neumann condition: 

^=0 x = o,b (15) 

3x 

Likewise the normal current j to the insulator spacing between electrodes must 
also vanish. Using equation (10) and noting from equation (6) that the velocity vanishes 
on the boundary yields the condition 

+ $£ = 0; z = o,h; 2 <y <a (16) 

3y 3z 2 

The final boundary condition relates the cell entrance region (y = 0) to the cell exit 
region (y = a), shown in figure 1, by a periodic condition, namely, 

<p(x, a, z) = <p(x,o, z) + <p E (17) 

where is the Hall potential. 


Transforming to Nondimensional Variables 

The Hall potential has ideal maximum that can be found from equation (9) when 
u = (u) , j =0 (open circuit), and dcp/dz = 0 (no local Faraday field). In this case, 
equation (9) becomes 


= ( u ) B/3 

3y 

Integrating over one cell length gives the ideal maximum: 

(<P H ) =<u>B/3a 

v 7 max 


(18) 


(19) 


The coordinates (x,y,z) are transformed with respect to the channel dimensions 
(b, a, h) (see fig. 1). Hence, 


x = b| 


( 20 ) 
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y = ai7 


( 21 ) 


z = H 


( 22 ) 


The ratio of channel height to width is given the symbol 


h 

V = — 
b 


(23) 


And the ratio of channel height to length is given the symbol 


h 

v = - 

a 


(24) 


The potential function <p is nondimensionalized by the substitution 


cp = (u) Bln// 


(25) 


The velocity function u is nondimensionalized with respect to (u) ; that is, 


u = (u>U 


(26) 


With these definitions and equation (6), the normalized function U becomes 


U = 


[?d - o] 1/m m - o] 1/m 


f 1 f 1 [1(1- %) ] 1 /m [r(i- 0 ] 1 /m d|<i t 
Jq Jq 


(27) 


Figure 2 is a plot of equation (25) at t = 1/2. Equation (11) becomes 

0£ 2 07J 2 0£ 2 ^ 


(28) 


The boundary conditions (eqs. (14), (15), (16), and (17)) become 


i// = 0; K =0,1; 0 =£ t/ < - 

2 


( 29 ) 
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0 ; | = 0,1 

dt 


(30) 




+ — =0; ? = 0,l;i<tj<l 
877 9? 2 




The maximum value of i//jj is then 


( ^h) “ ^ 

' 'max y 

where equations (19), (24), and (25) were used. 

The Hall current Ijj is nondimensionalized by letting 

Ir = o’ (u) BbhJjj 


(31) 

(32) 


(33) 


(34) 


then 


J H “ + ‘ 


1 + /3 


[ f (-y^+P^+l auWd? 
J 0 ^0 V dr l K / 


where the integral is computed at 77 = 3/4. 

Similarly, the Faraday current Ip is nondimensionalized by defining 

Ip = cr(u)BabJp 


(35) 


(36) 


then 


J F = 


--J -f 1 fU’-t 

1 + e 2 Jo Jo ' dT > 


77 + 


i^+ u)d£ d77 

9? / 


The integral, equation (35) is carried out over the plane £ = 1/2. 


(37) 
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Finite Difference Form of the Equations 


Each of the normalized coordinates (£ , t? C) spanning the region in figure 1 are 
divided into C - 1 equal line segments. 


K - 1 

(38) 

C - 1 


L - 1 

(39) 

C - 1 


M - 1 

(40) 


C - 1 


where (X,L,M) are integers ranging from 1 to C. The space is therefore subdivided 

q 

into a mesh having C points and the index of each point is given by (K, L, M) . 

The potential i//(i; , 17 , C) is replaced by an array A(K,L,M); and the increments 
Aij , A 77 , A£ are equal to 1/(C - 1). 

Three point Lagrangian differentiation formulas (see p. 914, ref. 7) are used to 
approximate the partial derivatives. 


A(K + 1 , L , M) + A(K - 1 , L , M) - 2A(K,L,M) 

3? 2 ( A |) 2 

A(K,L + 1 ,M) + A(K,L - 1,M) - 2A(K,L,M) 
d V 2 (At ?) 2 

9 2 i// M A(K,L,M + 1 ) + A(K, L, M - 1 ) - 2A(K,L,M) 
a ? 2 (A ?) 2 


(41) 


(42) 


(43) 


Substituting these approximations into the differential equation (eq. (28)), and solving 
for A(K,L,M), the potential at the point (K,L,M) can be written in terms of its neigh- 
boring points: 
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1 


A(K,L,M) = 


2 [l + r 2 + JLT 2 (1 + ^ 2 )] 


]i 2 (1 + 0 2 ) [A(K + 1, L, M) + A(K - 1, L, M)] 


+ y 2 [A( K , L + 1 , M) + A(K, L - 1 , M)] 


+ A(K, L, M + 1) + A(K, Xj, M - 1) + 


1 

(C - 1) 



( 44 ) 


There is, of course, no need to approximate BV/dC because this differentiation can 
be analytically found from equation (27) and does not contain the potential function. 

The boundary condition 3i^/3| = 0, (eq. (30)), which holds at 4 = 0 and ij = 1 (the 
sidewall insulators) must be handled in a way that insures that the approximation does 
not involve points beyond the region of solution (e.g. , points having negative or zero 
indices). At 4 = 0 (K = 1), 


3£ a -3A(K, L, M) + 4A(K + 1, L, M) - A(K + 2, L, M) = Q (45) 

34 2(A?) 

This gives the finite difference form of the boundary condition: 

A(K, L,M) =I[4A(K + 1,L,M) - A(K + 2,L,M)] , K = 1 (46) 

3 

At £ = 1 (K = C) the approximation is 


d± „ A(K - 2, L, M) - 4A(K - 1, L, M) + 3A(K, L, M) 
34 2(A£) 


from which 


A(K, L, M) = - 1 - [4A(K - 1, L, M) - A(K - 2, L, M)] , K = C (48) 

3 

Equation (45) is called a forward difference formula because it involves points forward 
of the point under consideration. While equation (47) is called a backward difference 


13 


formula. There is also a central difference formula: 


d± - A(K + 1,L,M) - A(K - 1, L, M) (4g) 

9? 2(A£) 

This form obviously does not involve the point A(K, L, M). In considering the boundary 
conditions at the electrode spacing insulators (eq. (31)) at C = 0 and tf = 1, a forward 
difference formula for dxp/dt must be used at f = 0 and a backward difference at 
if = 1. But equation (31) also involves the derivative dx}//dT} for which a central differ- 
ence formula could be used. However, it was found that the problem would often fail to 
converge unless a difference formula involving the point A(K,L,M) was used. The 
reason for this behavior is not clear but may be related to the type of boundary condition 
given by equation (31), which is neither a Dirichlet nor a Neumann type. Because of 
this behavior it was decided to use a forward difference at If = 0 for both dx(//d£ and 
dxjy/dri and a backward difference for both at £ = 1. Therefore, at if = 0 (M = 1), 

(1/2) < i) < 1, and [(C + l)/2 < L < C] equation (31) yields 

A(K, L, M) = 1 f/3y[4A(K, L + 1,M) - A(K, L + 2,M)] 

3(1 + /3y) L 

+ [4A(K, L, M + 1) - A(K, L, M + 2] | (50) 

It should be noticed that, at L = C - 1, equation (50) contains the point 
A(K, C + 1, M), which is beyond the range of the mesh points. Points such as these, 
that overrun the L range, are easily handled by using the periodicity condition. 

At if = 1 (M = C), (1/2) < 77 < 1, and [(C + l)/2 < L < C] 

A(K, L, M) = l f/3y[4A(K,L - 1,M) - A(K, L - 2,M)] 

3(1 + /3y) 1 

+ [4A(K, L, M - 1) - A(K, L, M - 2)]| (51) 

The remaining boundary conditions (eqs. (29) and (32)) are straightforward Dirichlet 
conditions and will not be rewritten here. 


Method of Solution 

Because the the large number of points involved in a three-dimensional problem 


14 



(for this study C = 17 and C 3 = 4913 = number of points), a relaxation method must be 
used to find the A(K, L, M). Direct methods such as a matrix inversion technique would 
require much more computer storage than is now available. On the other hand, relaxa- 
tion methods are notorious for oscillatory behavior. This problem can be handled by the 
dynamic relaxation method in which a damping factor is introduced. 

The method used in this report is a modification of the dynamic relaxation method 
described in reference 8. In that reference an artificial time -dependent problem is 
added. Equation (1) would then be modified to 

div 7 = D ^2 + (52) 

3t 3,2 


where D is a chosen damping factor. 

However, for this study, the time dependent portion will be introduced as an 
auxiliary equation that is solved simultaneously; that is, 


div j = 0 


(53) 


D 9£ + 9_£ _ q 
0t at 2 


(54) 


In terms of the normalized potential, the derivatives in equation (54) are written: 


01// a A i+1 (K, L, M) - A*~^(K, L,M) 


3t 


2(At) 


(55) 


0^i// a A i+1 (K, L, M) +A 1 ~ 1 (K,L,M) - 2A 1 (K,L,M) 


i-1 / 


dV 


(Aty 


(56) 


The superscripts refer to the iteration number. With these substitutions equa- 
tion (54) is solved for A 1+1 (K,L,M): 


A i+1 (K, L,M) = 1 


l + ^t) 

< 2 / 


2A i (K,L,M) - A i-1 (K,L,M) |l -5^ 


(57) 
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The i** 1 value of A(K,L,M) is the current value as calculated from equation (44) 
(which is, of course, the reduced form of eq. (53)). Only one complete set of values 
A(K,L,M) is stored. The set is continually overwritten by the latest calculations. This 
technique is called Liebman’s method (see, e.g. , p. 270 of ref. 9) and is the most 
economical use of computer storage for large iterative problems of the type considered 
in this report. 

At the beginning of each calculation for the potential at (K, L,M) the quantity V is 
assigned the present value of A(K,L,M), so that 

V s A i_1 (K,L,M) (58) 

A new value A X (K,L,M) is then calculated from equation (44), and immediately 
overwritten with A* + *(K,L,M) from equation (57). When the deviation 

|A 1 + 1 (K,L,M)-V| Jg 1 (59) 

|A i+1 (K,L,M)| 16 


for all (K,L,M) the problem is regarded as having converged. 

No extensive surveys of the range of the quantity DAt/2 has been made. However, 
it was found that with a (DAt/2) =0.1 convergence occurs at about 200 iterations. The 
problem calculated in reference 8 had 1064 points and requires 92 to 160 iterations de- 
pending on boundary conditions. 

Preloading a problem with the data from a previous case has a useful effect on the 
rapidity of convergence. Rarely does a preloaded case require more than 100 iterations 
and generally about 30 to 50. This is, of course, dependent on the choice of the initial 
data. 


RESULTS 

The equipotentials in various planes are plotted for the case of p = y = 1 (cubic 
geometry). These results are typical of the structure of the potential function. Two 
cases are presented: (1) Hall potential equal to 0. 25 (approximately midload line) and 
(2) short circuit (no Hall potential). In both cases a parabolic velocity profile is used 
(m = 1). 

Figure 3 is a plot of the equipotentials in the 4 = 1/2 planes for a Hall potential 
of 0. 25. One set of electrodes are located between L = 1 and L = 9, and a second 
electrode pair begins at L = 17. A region of high electric field is located at the trailing 
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Equipotential, 


4 > 



Figure 3. - Equipotentials in (L, M) plane at centerline (K * 9). 


edge of the lower electrode, the electron emitter, at L ^ 9. Another point of field con- 
centration is formed at the leading edge of the top Hall electrode, the electron collector, 
at L =* 16. These high field regions are characteristic of all the equipotentials that 
were found. The potential variation along the separating insulator is similar to that found 
in the two-dimensional studies of references 1 and 2. 

The points (5,9) and (13,9) have inversion symmetry properties in the (L,M) plane. 

If an orthogonal cartesian coordinate system is constructed with either point as origin, 
then a double reflection (first through one axis and a second reflection through the other 
axis) will leave the potentials unchanged. 

Figure 4 shows views of the equipotentials axially down the channel. The velocity 
is directed into the plane of these figures. Figure 4(a) is the entrance region of the cell 
(L = 1); part (b) midway across the electrodes (L = 5); part (c) at the trailing edge of 
the electrodes (L = 9); part (d) midway at the insulators (L = 13); and part (e) at the 
exit plane of the cell (L = 17). Note that the entrance and exit planes differ only by the 
Hall potential as is required by the boundary conditions. 

Figure 5 is the short-circuit counterpart of figure 3. And figure 6 shows the short- 
circuit counterparts of figure 4. The slight asymmetry that is shown in the short- 
circuit plots can be eliminated by tightening the convergence requirement (eq. (57)). 

This, of course, costs computer time and is only of aesthetic value because the integral 
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1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 

Distance parallel to magnetic field, K 

(b) Midway across electrode (L = 5). 

Figure 6. - Short circuit equipotentials in (K, Ml plane. Dimensionless 
Hall potential, 0. 












for the Hall current is affected only in the fourth significant figure if the convergence 
factor in equation (59) is changed to 1/32. 

Figures 4 and 6 are graphic displays of the three-dimensionality of the problem. 
Without a velocity profile, equation (28) could be reduced to a two-dimensional problem 
by separation of the £ variable. Because of the boundary condition 3 i/'/s 4 = 0 at the 
insulators, the potentials in the (£,?) or (K,M) plane would be straight lines. 

Figure 7 is a plot of Hall potential against Hall current with m (the exponent in the 
velocity function) as a parameter. Figure 8 is a plot of the Hall current for the cross 



Figure 8. - Integral of Hall component of current density as function of 
axial distance. Dimensionless Hall potential, 0. 25. 


sectional planes L = 1 to 17. It illustrates the reason for choosing the plane L = 13 
for the definition of the Hall current. Figure 9 is a similar curve of the Hall potential 
as a function of the absolute value of the Faraday current with m as a parameter. 
According to reference 5, these load lines can be independent of the velocity profile in 
the direction parallel to the magnetic field. The effect shown in figures 7 to 9 must be 
due to the velocity profile in the Faraday direction. 

The dependence of figures 7 to 9 on the velocity profile is contrary to the results 
reported in references 3 and 10. According to these references the load lines are inde- 
pendent of profile. Reference 3, however, posits an arbitrary current dividing line, 
through which no current can flow in the Hall direction. Reference 10 seems to rest on 
the assumption that the open-circuit Faraday voltage is defined by letting the Faraday 
current density vanish identically . This is not realistic because experimentally it is 
only the total integrated Faraday current that can be set equal to zero. 

Figures 7 to 9 can be compared with the two-dimensional analysis of references 1 
and 2. For the short-circuit case (zero load factor) the quantity X defined in refer- 
ence 2 is equal to the normalized Faraday current plotted in figure 9 of this report. At 
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zero Hall current, figures 7 and 8 show that the Hall potential is 0.47 to 0. 51. Figure 9 

shows that this is equivalent to a normalized Faraday current of 0.625 to 0.663. The 

value of X from reference 2 for the same geometry and Hall parameter is about 0.6. 

Reference 1 is in agreement with reference 2 for small Hall parameters (/3 — 1). 

o 

Figure 10 shows the variation of J along the axial distance L at the midpoint of 

the channel (K = 9). A Hall potential of 0.25 and a turbulent profile (m = 7) was used in 

o 

the computation. Regions of high J heating are found at the trailing edge of the cathode 
(lower electrode) and at the leading edge of the anode (upper electrode). Along the 
centerline of the channel (K = 9, M = 9) the ohmic heating is nearly constant with axial 
distance L. 

A cross plot of figure 10 in the Faraday direction M at L = 5, 9, and 13 is shown 

in figure 11. The station L = 5 is midelectrode; L = 9 is at the trailing edge of the 

2 

electrodes; and L = B is midinsulator. The plot shows that J is diminished at the 
channel walls (electrode -electrode spacing insulator) except for the hot spots near the 
trailing edge of the cathode. 

O 

Figures 12 to 14 show the variation of J at the same stations (L = 5,9,13) be- 
tween the sidewall insulators. These plots show that the walls are regions of greatly 
reduced ohmic heating. 

9 

Figure 15 shows the variation of J along the channel midplane (K = 9) for a para- 
bolic velocity profile (m = 1) . This figure is calculated at the same conditions as fig- 
ure 10. The only difference being the velocity profile. The same general characteris- 
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Axial distance, L 

Figure 10. - The square of normalized total current density J 2 distribu- 
tion as function of axial distance for turbulent profile (m « 7). Mid- 
channel (K = 9); dimensionless Hall potential, 0. 25. 



Figure 11. - Distribution of the sauare of normalized total current den- 
sity J 2 as function of Faraday distance for turbulent profile (m = 7). 
Midchannel ( K = 9); dimensionless Hall potential, 0. 25. 



)/(o(u> b) 2 j2 . (j2 + j2 + j|)/(a <u) b) 


Distance in 
Faraday direction. 



Figure 12. - Distribution of the square of normalized total current den- 
sity as function of distance across channel at midelectrode (L = 5) 
for turbulent profile (m = 7). Dimensionless Hall potential, 0. 25. 


Distance in Faraday direction, 
M 



Figure 13. - Distribution of the square of normalized total current den- 
sity J z as function of distance across channel at trailing edge of elec- 
trode (L = 9) for turbulent profile (m = 71. Dimensionless Hall potential 
0. 25. 
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TABLE I. - COMPARISON OF INTEGRATED THREE-DIMENSIONAL EFFICIENCY WITH LOCAL EFFICIENCY 


Velocity 

profile 

parameter, 

m 

Profile 

Dimen- 

sionless 

Hall 

potential , 

Dimen- 

sionless 

Hall 

current, 

J H 

^H J H 

Dimensionless 
ohmic heating, a 

- ft J 2 d£ drj d? 

v J 

Useful 

efficiency * 5 

Local 

efficiency 0 

Hall load 
factor, 

k h 

1 

Parabolic 

0.25 

0.188 

0.047 

0.599 

0.073 

0.163 

0.532 

7 

Turbulent 

.25 

.219 

.054 

.457 

.106 

.168 

.492 


a Ohmic heating = - J 3 2 dx dy dz s abhor ((u) B ) 2 y &K ■ 


Useful efficiency 


/T 


• E dx dy dz 


/ 


3 ■ E dx dy dz + 


a J 


dx dy dz 


-Ref. 11; eq. 14.17. 


^H J H 

^H J H + i f j2d? d?,dlr 

y J 


The useful power efficiency is 7.3 and 10.6 percent, and this is compared with the 
local efficiency of 16. 3 to 16. 8 percent (ref. 11). Thus, volume effects reduce the effi- 
ciency by about a factor of two. 

The last column in table I gives the integrated ohmic heating for both velocity pro- 
files. The integration is scaled by a factor abhcr((u) B) 2 . Multiplying by this factor 
will give the ohmic heating per unit volume. 


CONCLUSIONS 

Electron heating of a plasma moving through a magnetic field is directly propor- 

o 

tional to the square of the current density (J ) . Preferential electron heating causes a 
nonequilibrium conductivity in the plasma. 

The phenomenon called boundary -layer shorting presumed that the channel walls 
were regions of such enhanced conductivity. This view was held because the reduced 
velocity profile near a wall would result in reduced local electric fields. These local 
fields drive the current against the induced electrode potentials. Near the walls, how- 
ever, the local fields are smaller, and it was expected that a large current would flow 
in the opposite direction and short the adjacent electrodes. 
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As a result of solving for the three-dimensional potential and current distributions 
in an MHD generator with various specified velocity profiles and constant conductivity 
the following conclusions are made: 

1 . The sidewall insulators of the channel are always regions of lesser heating 
than the central portion of the channel (see figs. 12 to 15). Nonequilibrium conductivity 
is therefore more likely in the central portion than along the sidewalls . 

2. The electrode wall insulators (segmentation insulators) are also regions of lesser 
heating than in the central portion of the channel, except at the trailing edge of the 

cathode and at the leading edge of the anode. The gradient of these regions is so steep 

9 

that, slightly further along the insulator, away from the electrode, the J heating is 
insignificant (see figs. 10 and 15). Unless these hot spots produce "pockets" of highly 
conducting plasma that migrate along the insulator to the adjacent electrode, it is diffi- 
cult to see how shorting can take place. 

3. Different velocity profiles provide differently sized regions near the walls in 

which the local fields are reduced and in which, presumably, the shorting phenomenon 

2 

might be observed. For all the profiles used in this study, the calculated J distribu- 
tions are substantially the same in that the walls always have less heating. 

4. Surface effects on the insulator material, such as coating with seed material, 
are outside the scope of this study. Such effects could explain observed shorting and 
still be consistant with the results of this study. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 24, 1971, 
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